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ABSTRACT 

We report the results of a joint Chandra-HST study of the X-ray binary pop- 
ulation in the massive, high-density globular cluster NGC 6388. NGC 6388 has 
one of the highest predicted X-ray binary production rate of any Galactic clus- 
ter. We detected a large population of 61 Chandra sources within the half-mass 
radius with hx > 5 x 10^° erg s^^ From the X-ray colors, luminosities, (lack of) 
variability, and spectral fitting, we identify five as likely quiescent low-mass X-ray 
binaries. Due to the extremely crowded nature of the core of NGC 6388, finding 
optical identifications to Chandra sources is challenging. We have identified four 
blue, optically variable counterparts to spectrally hard X-ray sources, evidence 
that these are bright cataclysmic variables (CVs). One showed variabihty of 2 
magnitudes in V, indicative of a dwarf nova eruption. One other likely CV is 
identified by its X-ray spectrum (partial covering with high Nh) and strong vari- 
ability, making five likely CVs identified in this cluster. The relatively bright 
optical magnitudes of these sources put them in the same class as CVl in M15 
and the brightest CVs in 47 Tuc. 

Subject headings: globular clusters: individual (NGC 6388) — X-rays: binaries — 
novae, cataclysmic variables 
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1. Background 



Studying the relationship between X-ray binaries (XRBs) and globular clusters 
continues to provide interesting clues into cluster dynamical evolution. These clusters 
have long been thought to be major producers of XRBs, since dense cluster c ores promot e 



the type of dynarn ical interactions that lead to tight binaries such as XRBs (j Clark 



Pooley et al. 



1975 



20031 ) . These interactions include exchange interactions where a massive stellar 



remnant such as a heavy white dwarf or neutron sta r has a close encounter with an existing 



binary, replacing one of the other stars (Hills 



a main-sequence star (IFabian et al 



19761 ). tidal capture by a heavy remnant of 



19751 ). or collisions of neutron stars with red giants. 



eventually prod ucing ultracompact binaries with neutron stars accreting from white dwarfs 



( Verbunt 



19871 ). Cataclysmic variables (CV s) and low-mass X-ray bin ary (LMXB) systems 



are common products of these interactions (jivanova et al. 



2006 



20081 ). In globular clusters. 



a significant numbe r of millisecond pulsars (MSPs) have been identi fied, such as in 47 Tuc 



(ICamilo et al 



20001 ). Such systems are the likely progeny of LXMBs fArchibald et al. 



20091 ). 



Since mass segregation causes the more massive populations of binary stars and degenerate 
remnants (for instance, MSPs) to be more centrally concentrated, these objects will be 
strongly tied to the dynamical evolution of the cluster. Most significantly, this increases the 
likelihood of further interactions of binaries in the core which ca n then halt th e progress of 



dynamical relaxation, supporting the core against core collapse (jFregeau 



20081 ). 



The known distances, reddenings, and ages of globular clusters allow for an analysis of 
populations of CVs and LMXBs that is not possible for field objects. Using the Hubble 
Space Telescope (HST) in conjunction with the Chandra X-ray Observatory (Chandra) 
large numbers of X-ray sour ces have been identified in clusters and matched to optical 



counterparts. For instance. 



Heinke et al. 



( l2005l ) detected 300 X-ray sources down to 



Lx ~ 8 X 10^^ erg s ^ within the half-mass radius of 47 Tuc, 105 of which could be 
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securely classified as quiescent LMXBs (qLMXBs), CVs, MSPs, or chromospherically active 
binaries (ABs), mostly by using HST to identify and characterize optical counterparts 



( lEdmonds et al 



2003a 



b|). Similar identifications were also made of 59 of the 79 X-ray 



sources detected (with limi ting Lx ^ 1 0^^ er g s ^) in NGC 6397, a cluster that has 



undergone core-collapse, by 



Cohn et al 



(120 lOl ) in a joint HST- Chandra study. 



Hig her numbers of int eractions in cluster cores lead to larger populations of X-ray 



sources fjPooley et al 



2003h. This interaction rate can be simply estira ated by the encounter 



rate, F oc pIt^/vq (iVerbunt fc Hut 



1987 



Maccarone fc Peacock 



20111 ). This quantity is a 



measure of how often stellar interactions will occur as a function of the cluster's central 
density, core radius and central velocity dispersion. The numbers of X-ray sources in 
a cluster, particularly the qLMXBs, MSPs, and brighter CVs, vary directly with this 



encounter rate, implying strongly that 



these populations in globular clusters (jPooley et al. 



dynamical interactions are the 



2003 



Heinke et al 



20061 ). while fainter ABs seem to be largely primordial in origin (IBassa et al. 



jrimary source of 



2003 



Poolev fc Hut 



2004h . 



1.1. NGC 6388 

NGC 6388 is a particularly good candidate for the study of the dynamical origin of 
XRBs in dense-core globular clusters due to its larg e central den sity (10^'^^ LqPc"'^) and 



core radius (7.2" = 0.35 pc at 9.9 pc from the Sun) (IHarris 



1996 



These give it one of the 



largest values of P of any globular cluster in the Galaxy. With such a large encounter rate, 
we expect it to host a rich populat i on of LMXBs and CVs, and thus to find a large number 



of X-ray sources. In 



Lugger et al 



(119871 ). the U-band surface brightness profile, which is 



dominated by horizontal branch stars, was found to be well fit by a single mass King model. 



http://physwww.physics.mcmaster.ca/~harris/mwgc.dat for 2010 revision. 
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suggesting that the core is hkely being supported by a large number of primordial binaries. 
It is possible that some of these primordial binaries could be detected now as ABs. 

NGC 6388 also has several features beyond its large encounter rate that make it 
a cluster worth investigating. For inst ance, for a cluster of its metallicity (as high as 



[Fe/H] 



0.7 flWallerstein et al. 



20071 ). it has a rather atypical extreme blue horizontal 



branch (EBHB) in additi on to a blue horiz ontal branch (BHB), and a very populated red 



horizontal branch (RHB). 



Rich et al. 



(119971 ) investigated whether this could be the result of 
dynamical interactions in the core enhancing the BHB via stripping of the envelope mass. 
However they found no difference in the radial distribution of the RHB and BHB stars 
and took this as evidence against the dynamical argument. Also intrig uing are the blue 



hook st ars (hotter and less lumino us than EBHB stars), investigated in 



Dalessandro et al. 



(120081) 



Noyola fc Gebhardt 



brightness profile and 



(120061) report e vidence for a modest central cusp in the surface 



Lanzoni et al 



( 120071 ) see a modest central surface density cusp 



in the inner arcsec, which are interpreted by those groups t o be eviden c e for a central 
intermediate-mass black hole (IMBH) of mass 5.7 x lO^M©. 



CsehetaL 



mm show that 



their radio nondetection gives a n upper limit to the mass of a c entral IMBH that i s less 



than half the mass predicted by 



Lanzoni et al. 



(120071 ). Recently 



Liitzgendorf et al. 



used s tellar kinematics to argue for a central IMBH with a mass of 17 ±9 x IO^Mq. 



(1201 ih 



Miocchi 



(j2007l ) suggests that tidal stripping of red giant envelopes by this putative IMBH may be 



producing the significant population of EBHB stars. 



VerbuntI (120011 ) analyzed two ROSAT PSPC observations of NGC 6388 from 1991 
and 1992, identifying an X-ray source at Lx (0.5-2.5 keV)=6 x 10^^ ergs/s just outside the 
half-mass radius. We do not detect a source in the smaller of Verbunt's two error circles. 
This could indicate a transient was active during the ROSAT observations, and faint in 
quiescence during the Chandra observations. It also could indicate that Verbunt's ROSAT 
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astrometry was incorrect, since the luminosity is close to our measured total Lx for the 
cluster. Less sophisticated X-ray analyses of our Chandra observation and a 2003 XMM 
observation, with the goal of constrainin g the presence of a n intermediat e - mass black hole 



at the cluster center, were performed by 



Cseh et al. 



(12010h : lNucita et all (12008[ l. 



Evidence that NGC 6388 contains at least one transient X-ray binary (in agreement 
with predictions of its high stellar encounter rate) was provided by IGR J17361-4441, a 
hard X-ray transient located in NGC 6388. It was detected by INTEGRAL on August 11, 
2011, with follow up observations by Swift, RXTE, XMM, ATCA, and Chandra. It had 
a peak output in t he 0.5-10 keV energy band of Lx ~ 5 x 10^^ ergs/s (IBozzo et al 



Nucita et al. 



2011 



2OI2I ). A short Chandra target of opportunity observation provided a precise 



position, inconsistent with any of the cluster X-ray sources we pre sent below, a " 
upper limit on the transient's quiescent luminosity of 10'^^ ergs/s (jPooley et al 



lowing an 



201lh . In 



the last Swift observation of Nov. 4, 2011 (before solar proximity forbade observations), this 
transient was still active at similar luminosities, but the X-ray em ission from the clu ster 



was consistent with the normal quiescent output by Jan. 29, 2012 fiBozzo et al. 



2OI2I). 



2. Data 



2.1. Chandra Observations 



We obtained a 45 ks ACIS-gf| observation of NGC 6388 on 2005 April 21/22, with the 
center of the cluster placed on the S3 chip. We chose the S3 chip for its high sensitivity and 
spectral resolution. We used the timed exposure mode wit h the very fa int telemetry format. 
The moderate reddening of NGC 6388 {E{B -V) = 0.37, 



Harris 



19961 (2010 version)) gives 



a moderate hydrogen column density, Nh ~ 2.2 x 10^^ cm ^. A detection limit of 2 counts. 



^Advanced CCD Imaging Spectrometer /S- Array 
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for the 45 ks ACIS-S exposure, corresponds to a flux of /x ~ 4 x 10 erg s ^ cm ^ in tlie 
0.5-6.0 keV energy ba nd for a pow er-law spectrum with a photon index of 2. Assuming 



a distance of 9.9 kpc ([Harris 



1996 



(2010 version)), this limiting flux corresponds to a 
luminosity of Lx ~ 5 x lO^'' erg s~^. We are complete to ~ 5 counts, except in the core 
and near bright sources. We are likely complete to ~ 10 counts except very near to bright 
sources. In Fig. [T] we plot the location of the identified X-ray sources on the X-ray image. 
We constructed a false-color representation of the sources. Fig. [2l where the variety of soft 



and hard X-ray sources is easily noted. 



2.2. HST Observations 

Our optical analysis is based on ACS/HRcfl imaging obtained in 2003 and 2006 
through the F555W (hereafter I/555) and F814W (hereafter Jsm) filters (GO-9835 and 
GO- 10474, PI: Drukier) and imaging obtained in 2006 through the F330W (hereafter 
[/330) filter (GO-10350, PI: Cohn). We refer to the filters used in a way that suggests 
similarity to the Johnson system, but with subscripts to indicate that the obtained data is 
not transformed to that system, especially true in the F330W filter, which is both bluer 
and narrower than the Johnson U filter. Table [T] lists the data sets, dates, filters, and 
exposure times. The observations used heavy dithering, with integer pixel plus subpixel 
offset, in order to allow for bad pixels and to increase the effective resolution through 
drizzle reconstruction of the frame. Rather than using the standard STSDAS multidrizzle 
algorithm, we used images created through a similar procedure, applied by J. Anderson. 
The resulting drizzle-reconstructed [/330 image is shown in Fig. |3l We note that while the 
dithering substantially improves the effective resolution by about a factor of two, the core 

■^Advanced Camera for Surveys/High Resolution Channel 
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region of NGC 6388 remains very crowded, with a considerable amount of psf overlap. 



3. X-ray Analysis and Results 

We reprocessed the Chandra data using CIAO 3.2, CALDB v. 3. 0.2, and standard 
procedureqj. This involves identifying and excluding bad pixels, utilizing information kept 
in the VFAINT telemetry format to reduce background, using accurate time-dependent 
gain and CTI corrections, removing pixel randomization, and selecting for grade and status. 
No background flares were seen, so we used the full 45.16 ks of live observing time. 



3.1. X-ray source identification 



We used the CIAO wavdetect program on Chandra images of the cluster within the 
half-mass radius, in the 0.3-6 and 0.3-1.2 keV energy ranges, to detect 61 sources. The 
source positions were refined, and net counts and spectra extracted, using the IDL package 
acis_extract (Broos et al. 2010); see Table 3. Using the net counts, and a conversion from 
counts to flux for each of 7 narrow bands, we computed an X-ray color- luminosity diagram 
(Fig. H]). Five sources (marked as squares) have relatively soft X-ray colors while being 
relatively luminous {Lx > 10^^ erg s^^). We plot the expected positions in this diagram 
of a 10-km radius, 1.4 MQueutron star (NS) with low B-fleld and a hydrogen atmosphere 
(NSATMOS, Heinke et al. 2006). The flve soft sources are consistent with this line, or 
slightly harder (which may indicate a spectrum containing both a NS atmosphere and 
a harder component). By comparison w ith similar diagrams in other globular clusters 



(IHeinke et al 



2003 



Poolev fc Hut 



20061 ). we infer that the harder sources (consistent with 



^http:/ /cxc. harvard.edu/ciao/threads/all. html 
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e.g. a power-law of photon index 1.5) are predominantly CVs. Below we identify evidence 
for the CV nature of five of them (labeled with blue triangles). 



3.2. X-ray Spectral Analysis 



Below we discuss X-ray spectral analysis of those sources with more than 80 counts. 
Our principal objective was to identify which sources have spectra typical of quiescent 
LMXBs, vs. which sources have spectra more commonly associated with cataclysmic 
variables. We used spectral fits o f a power-law, MEK AL (hot thermal plasma, with the 



cluster abundance of -0.6, or 25%; 



Liedahl et al. 



19951 ) which is of ten appropriate for fitting 



Heinke et al. 



20061). We 



CVs, and NSATMOS + power-law (typical of quiescent LMXBs, 
performed spectral fits in which we binned all spectra, using 20 counts/bin for brighter 
sources and 10 cts/bin for sources below 150 counts, and directly on the unbinned data 
using the C statistic. The results from these fits were similar, but the unbinned spectral 
fitting allowed better constraints, particularly on the fraction of flux in the (above ~2 keV) 
power-law component in NS+power-law spectral fits. Results from the unbinned spectral 
fits are presented in Table HI while example binned spectral fits are presented in Fig. [5l For 



all fits we imposed a minimum Nh of 2.2 x 10 cm , the value inferred from the reddening. 



3.2.1. Strong candidates for quiescent LMXBs 



Five sources (CXI, CX2, CX3, CX8, and CXIO; Table Hj) require very large photon 
indices (> 3) when fit by a power-law model, or very low MEKAL temperatures (< 1.5 
keV) for that model. Spectral fits to NSATMOS models plus power-law components give 
excellent fits without need of the power-law component, whic h can be constrained to <30% 
of the 0.5-10 keV flux (for assumed photon index of 1.5, e.g. 



Cackett et al. 



2010). In four 
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of them, the power-law component can be constrained to < 13% of the flux. We take 
these spectral fits as strong evidence for a quiescent LMXB interpretation. Combined with 
the location of all five within or close to the cluster's core, and the lack of bright optical 
counterparts (which argue against foreground coronal sources), we conclude that these five 
are quiescent LMXBs. 

3.2.2. Hard X-ray sources 

Most of the remaining brighter X-ray sources have spectra consistent with typical 
CVs in other clusters-consistent with high-temperature (~5-10 keV) MEKAL models with 
absorption similar to the cluster value (Table Hj). However, CX4 has a very hard spectrum 
(photon index T = 0.9^q 4). This suggests a CV nature, as CVs tend to have hard intrinsic 
spectra, and have more intrinsic absorption than active binaries or millisecond pulsars. 
This source cannot be well fit by MEKAL models using standard absorption (see Table HI 
Figure |5]). It can be well fit by a simple absorbed power law, but with a photon index of 
0.301 

o.tly unusually flat for an intrinsic X-ray spectrum. A partial covering absorber model, 
commonly used for intermediate polar CVs, provides a good (59% "goodness") fit; using 
Nh (PHABS) of 2.2 X 10^^ plus a PCFABS component covering 92% of the emission with 
3.9 X 10^2 ^^-2 rjj^jg gj^gg j^rp ^ 3.7^^5 keV, and Lx = 2tJ:^ x lO^^ erg s'^ The high 
luminosity and hard spectrum suggestive of partial covering suggest that this source may 
be an intermediate polar CV. 

3.3. X-ray source variability 



We used the glvary program in CI AO flGregory &: Loreddll992l ) to search for variability 
in all X-ray sources with more than 30 counts. Only one source produced a significant 
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detection of variability, CX4, which produced an odds ratio of 10^^'^ and probabihty of 1.0 
of variabihty. Its hghtcurve (along with that of a typical nonvariable source, CX2, with odds 
ratio 10~^'^ and probability 0.043) is shown in Fig. [HI We note that the quiescent LMXBs 



in NGC 6388 do not show obvious variability, as expected i : 



reradiation of heat stored in the deep crust ( iBrown et al 



thei r X-ray emission arises fro m 



1998 



Ushomirsky &: Rutledge 



20011). Some quiescent 



(IFridriksson et al 



2010 



MXBs have shown clear variability in their thermal emission 



Cackett et al. 



20101 ). but the subset of quiescent LMXBs without 



a strong power-law spectral component have not yet shown evidence for such variation 



teeinke et al. 



20061) 



3.4. X-ray Source Distribution 

By examining the spatial distribution of the X-ray sources in NGC 6388, we may 
infer the masses of the objects fro m which the X-ray fluxes originate (see for instance 



the discussion in 



Cohn et al. 



20101). We compare the distribution of the X-ray sources to 
the cluster parameters derived from the distribution of stars used in finding the center of 
the cluster (see §4.31) . We expect the X-ray sources to be more centrally concentrated. 



and thus more m assive, assumirig the 



previous studies (iGrindlay et al. 



2002 



populations 



Cohn et al. 



are in thermal equilibrium. As with 



20101) we fit the different populations 



with generalized King models: 



S{r) = So 



1 + 



a/2 



where tq is related to the core radius Tc by 
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(2) 



Using the stars with V555 < 20.0, we find the optical distribution of NGC 6388 to be well 
fit by = 7.2 arcsec and a = —2. If we assume thermal equilibrium, the slope parameter 
for the X-ray sources ax is related to the slope parameter for the optical distribution by 



ax = 1 + q{ao - 1] 



(3) 



where q = Mx/M^ is the ratio of the characteristic mass of the X-ray sources to the 
mass of the objects that dominate the optical light distribution. We fit the King model 
using a maximum-likelihood algorithm. We then applied a bootstrap method, where 
1000 resampled populations are created by resampling from the original population until 
the total number is reached, replacing the chosen source back into the pool before each 
choice. We used the set of these resampled populations to determine uncertainties of the 
best-fit parameter values by fitting the King model to each resampled population and then 
analyzing the resulting distribution of the best-fit parameter values. The uncertainty in 
each parameter was set equal to half of the 68% spread of the best-fit parameter values 
about the median. For the fits we used the center of the cluster as determined from the 
optical data (see §4.3p . 

Since the cluster lies in a region where the Galactic ri dge emission is weak (1= 345.56° 



and b=-6.74° 



Harris 



1996 



Revnivtsev et al. 



20061). the 



20011) to 



(2010 edition); see the map in 
background contamination was only determined using Eqn. 1 in ( iGiacconi et al. 
estimate extragalactic contamination. We find that over the area covered by all the detected 
sources, there are likely 2-3 background source detections. In an effort to verify this value. 



we performed a n additional fit inc 



Eqn. 13. 4[ as in 



Eisner et al 



uding the additional parameter Bq for the background to 



(|2008[ ). This gave a result for Bq that was consistent with our 



previously calculated value. 
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The fit for the distribution of the X-ray sources by a King model profile along with a 
term for background sources is seen in Fig. [71 It is plotted alongside the distribution for 
the stars with V555 < 20.0 for comparison. From the fits (see Table [H]), we find the ratio of 
X-ray source mass, Mx, to the typical mass of stars with V555 < 20.0 to be g = 1.51 ± 0.12. 

We can infer the mass of t he stars used t o dete rmine the cluster parameters by 



comparison with 47 Tuc, which 



Catelan et al 



composition and age to NGC 6388. 



Heinke et al 



( 20061) h ave found to be very similar in 



(|2005[ ) found that for the range of stars in 



47 Tuc similar to those we used to determine the optical profile for NGC 6388, the average 
mass was approximately 0.88±0.05Mq. We use the approximate value of 0.9 Mq, to infer 
a mean X-ray source mass of 1.36 ± O.IIMq. This is consistent with a mixture o f mass ive 



Cohn et al. 



(2010) 



Heinke et al 



cataclysmic variables and low mass X-ray binaries. For comparison, 
report ed a mass of 1.14±O.14M0 for CVs in NGC 6397. And for 47 Tuc, 
( 120051 ) found a mass of 1.43±O.17M0 for all of the detected X-ray sources, and 1.31±O.22M0 
for the identified CVs. 

We also fit the identified qLMXBs, but ignored background contamination due to 
the small sample. This gave a value of g = 1.63 ± 0.31, which corresponds to a mass of 
approximately 1.47 ± O.28M0. While the difference in mass is not statistically significant, 
it is suggestive that the qLMXBs are more massive than the other X-ray sources (see Fig. 
E]). We also can compare the bright X-ray sources, defined as those with counts > 40, with 
the faint X-ray sources. The fits to these two populations can be found in Table |6] and their 
distributions are plotted in Fig. O 

We use a Kolmogorov-Smirnov (K-S) test on two samples to determine the likelihood 
of their radial distribution being drawn from the same population of stars. We find they 
are consistent with being from the same distribution at a probability of only 9%. While 
this does not represent a high level of statistical significance, it suggests that bright X-ray 
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sources may be more cen trally concentrate d than fainter ones, similar to that seen in NGC 



6397 for the bright CVs (ICohn et al. 



20101 ) ■ This may indicate that dynamical interactions 



scatter hard binaries into elongated orbits where they dim as they age. 



4. HST Data Analysis and Results 
4.1. Astrometry 

We put the ACS/HRC frames on the ICRS astrometric grid by first obtaining an 
astrometric solution for a ground-based frame, then using this frame to define secondary 
astrometric standards, and finally using these secondary standards to obtain a solution 
for the ACS frames. The ground-based frame, from the ESO archive, was taken with 
the Superb Seeing Imager 2 (SUS12) on the ESO New Technology Telescope (NTT) on 
2003 May 05 with a field of view of 5.46' x 5.75'. The primary astrometric standards are 
from the second US Naval Obse rvatory CCD Astrogr aph Catalog (UCAC2), which has 



an RMS uncertainty of 70 mas (jZacharias et al. 



20041 ). Consistent with this uncertainty. 



the astrometric solution for the ESO frame using 108 primary standards had an RMS 
uncertainty of 72 mas. Since the internal error of the ACS astrometric solution using 15 
secondary standards is about 13 mas, the overall accuracy of the astrometric solution of 
the ACS images is 73 mas. We took the corrected HST coordinates to be our fundamental 
coordinate system and applied a boresight correction to the Chandra source positions as 
described below in section 14.41 



4.2. The Color-Magnitude Diagram 



We used DAOPHOT-ALLSTAR to carry out PSF-fitting photometry on the two stacked 
V555 frames (2003 and 2006 epochs) and on the single stacked t/330 frame (2006 epoch). The 
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PSF was constructed using about 100 stars with the least crowding in the frames, with an 
iterative approach to neighbor subtraction. We constructed the PSF by first solving for it 
with a total radius that was relatively small: about twice the FWHM of the stars in the 
image. After each solution for the psf, a subtracted frame was constructed to clean the 
region around the chosen psf stars. A new psf with a larger radius was determined from this 
cleaner frame. This process was iterated until the change in the error to the psf fit became 
negligible. A total of 28612 objects were jointly detected in the 2006 V555 and U330 data. 

The magnitudes were calibrated to the STMAG system using the method described by 



Sirianni et al. 



(|2005[ ). Since the stacked images used for our PSF analysis had an arbitrary 
normalization, we constructed a second set of stacked images using multidrizzle. We then 
obtained aperture photometry from a subset of our stars from these images to the STMAG 
system. Specifically, we used the psf stars we had chosen for the DAOPHOT-ALLSTAR 
analysis. These magnitudes were converted to the STMAG system and then an offset was 
found between these derived magnitudes and the previously measured psf photometry. This 
offset was then applied to all the psf magnitudes in each of the filters and epochs. 

Figure [TO] shows the resulting (V555, f/330 — V555) CMD for the 2006 epoch - the epoch 
for which both U and V band imaging are available. The previously noted features of the 
CMD, discussed in §1.11 are apparent: the sloped red horizontal branch, the BHB, and the 
EBHB. Also we see a large blue straggler (BS) population, although the smearing out of the 
fiducial sequences makes it difficult to select blue stragglers near the main sequence turn-off 
(MSTO) with certainty. The MSTO lies about 3.7 magnitudes below the mean magnitude 
of the RHB and the MS extends about another 1.5 mag below the MSTO. Of note is the 
extension of the EBHB to about 1 mag below the MSTO and the presence of stars that lie 
below the MSTO and which are much bluer than the MS. As we discuss in §4.51 below, some 
of these stars appear to be counterparts to Chandra sources. It is useful to compare these 
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results to those of 



Catelan et al 



(120061 ). who have presented a deep HST-hased CMD of 
NGC 6388. They used data from the HST SNAP-9821 and GO-9835 datasets; the former 
was obtained with the ACS/WFC while the latter, which was obtained with ACS/HRC for 
the present program, provides higher angular resolution. Their analysis of these data with 
DAOPHOT-ALLFRAME produced a {V,B-I) CMD. While their overall CMD is relatively 
similar to ours, taking into account the different color indices used, we note that the EBHB 
in our CMD extends to an additional 2.5 magnitudes below the RHB, i.e. a total of 4.4 mag 
below the RHB. While the tip of the EBHB region is relatively sparsely popul ated, it raises 
many questions regarding stellar evolution (see, for instance, the discussion by 



Brown et al. 



2010|). 



4.3. Optical Radial Distribution 

The high resolution of the cluster core provides an opportunity to measure the position 
of the cluster center and analyze the radial distributions of the different populations of 
stars seen in the optical. To avoid problems with completeness, we use only stars brighter 
than 20 in V555 to determine the optical center of the cluster in the ICRS reference frame. 
This sample is used specifically to mirror the evaluation of the cluster center as described in 



Lanzoni et al 



(l2007f ). The centering method involved iterative centroiding over a circular 
aperture of radius 9". Using this sample of stars we find a center of RA = 17'^36'"17.185* 
± 0.12" dec = -44°44'06.44" ± 0.12". To find errors for the center, we use the bootstrap 
method as described in §3.41 We find that due to the small field of view for the ACS/HRC, 
the center cannot be determined to better than about 0.12 arcsec. This is because the 
cluster's core fills much of the field of view, giving little leverage for determining the center. 

In Fig. [HI we compare our computed center with other recently reported cluster 
centers. We find that our center differs by about 2-3 a from the centers determined by 
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Lanzoni et al.l (|2007| ) and iLiitzgendorf et al.l (1201 if ), lying somewhat to the north west. 
With our higher resolution, we were able to count more sta rs than either study, which may 



Noyola fc Gebhardt 



fl2006h was 



account for the differences. Note that the center derived by 
from a frame that was not corrected to a standard world coordinate system, which (as they 
point out) can be significant in HST fields. 

The stars used for finding the cluster center were then used to solve for the cluster 
parameters using maximum-likelihood routines to fit a generalized King model to the 
distribution (see Eqns. 1 and 2). Bootstrap statistics are used to determine errors on the 
fit. These stars give a value of 7.1 ± 0.2 arcsec for the core radius and —2.04 ± 0.08 for the 



slope of the 



respectively (IHarris 



all off. The se are consistent with previous solutions of 7.2 arcsec and -2.0, 



19961 ). which we use for comparison with other populations of interest. 



We fit several different stellar sub-populations in the same way as we fit the X-ray 
distributions (Table [6]). We see that the RHB core radius fit is consistent with the 
parameters previously found for this cluster of 7.2 arcsec. We also find it is has an 
a = —1.92 ± 0.33, which is consistent with the parameters for a relaxed system. A K-S test 
shows that the RHB distribution is consistent with the distribution of the other stars with 
V555 < 20.0 (Table [5]). Since HB stars will typically be less massive than the RGB stars 
that make up a majority of the the stars with V555 < 20.0, we assume they have not yet 
re-thermalized at their new masses. 

Fitting the blue stragglers (BSs) in the cluster shows them to be significantly more 
centrally concentrated than the stars with V555 < 20.0. We infer a mass of 1.80±0.27 times 
the mass of the stars with V555 < 20.0, giving a mass of approximately 1.62Mq for the BSs. 
When we plot these distributions, we can see how the BSs stand out from the stars with 
V555 < 20.0, while the RHB is consistent with that population (Fig. [T2l) . When we compare 
the BSs to the stars V555 < 20.0 with a K-S test, we find they are not the same, with a 
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probability of less than 1% of being drawn from the same underlying population (see Table 
E]). This is in agreement with the significant difference in the masses we find for the two 
populations. 

We also compare the distributions of the RHB, BHB, and EBHB stars to see if any 
one is more centrally concentrated than the others as a test for the influence of dynamical 
effects on the differences between these populations. K-S tests show that there there is no 
significant difference in the radial distributions of the RHB, BHB, and EBHB. Thus we find 



no difference in the masses of these thr ee groups o. 



previous findings for these populations (IRich et al. 



stars. This result is in agreem ent with 



1997 



Dalessandro et al. 



20081 ). 



4.4. Boresight Correction 

From the CMD, a very blue star was easily picked out. This star coincides with one of 
the brighter X-ray sources (source CX5), and provided a preliminary boresight correction 
to the X-ray positions. After applying this correction, we attempted to discover additional 
counterparts from their position in the color-magnitude diagram and proximity to X-ray 
source positions. But of the stars for which f/330 and V555 magnitudes were obtained 
through the standard procedures, there were none that could be identified in this way. 
Fortunately, three additional optical counterparts were identified by blinking the images in 
two bands (f/330 from 2006 against V555 from 2003 and 2006) and searching for blue stars 
near the X-ray positions. Using these four counterparts, a final correction of 0.105 arcsec 
in RA and -0.031 arcsec in Dec was applied to the Chandra positions, well short of the 
90% pointing error circle of ~ 0.6 arcsec. Once this correction was applied, we searched 
for additional counterparts. In the region near source CXll, a blue object was identified 
visually. Unfortunately it lies in the Airy ring of a much brighter star, preventing accurate 
photometry. We also extended our search to objects with a red excess in the f/330, V555, and 
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/840 bands in search of millisecond pulsars and active binaries. These additional searches 
yielded no other detectable optical counterparts. 

The small number of detected counterparts is not surprising considering that the 
factors that predict a large number of X-ray sources in NGC 6388 also make doing deep 
photometry challenging. Since the likely counterparts to X-ray sources are most likely to 
be centrally concentrated, the optical images were centered on the cluster. Because of the 
small field of view of the HRC on the ACS, the moderately large core radius represents a 
significant fraction of the entire field. As the core is also very dense, most of the optical 
counterparts must be found against a bright background of the core stars, making sky 
subtraction unreliable for objects more than about 3 mag below the MSTO. This also 
results in X-ray sources somewhat outside the core (such as CX4) falling outside of the field 
of view in the optical images. 

4.5. Analysis of Optical Counterparts 

We performed photometry for each of the identified optical counterparts using 
DAOPHOT-ALLSTAR on combined, drizzled frames for each of the epochs (Table [2]). One 
can see from the CMD (Fig. [TOl) that these four stars lie well to the blue of the main 
sequence (MS) and near the limit of detection in the V555 frames. For three of these 
stars, a light curve was constructed using aperture photometry performed on the original, 
undrizzled frames. The points were then further binned to provide a statistical data set 
with which we could test the variability of the optical sources. For the dimmest counterpart 
in the densest part of the field (source CX7), measuring in undrizzled frames was not 
possible in either the 2003 or 2006 epoch. In this case, we drizzled and combined the data 
into three images. For reference, finding charts are shown in Fig. [131 
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4-5.1. Source CX5 



CX5's counterpart is extremely blue (f/330 — V555=-1.2). The Airy disk from a brighter 
nearby star makes /814 measurements impossible. Our V555 light curves for 2003 and 2006 
(Fig. [T^ show variability of 0.6 and 0.8 mags (respectively), on timescales of hours. The 



2003 variability suggests sinusoidal variations with P> 8 h, but the amp 
high for ellipsoidal variability (typically 0.1 mags, e.g. Edmonds et al. 



itudeis extremely 



2mm- The 



companion's average brightness drops by 0.7 mags from 2003 to 2006. The si ze of these 



changes may suggest dwarf nova eruptions, which are typically 2-5 mags (e.g. 



Shara et al 



19961). although disk flickering cannot be excluded. Using the distance of 9.9 kpc from 



Harris! ( 1l996l . 2010 version) and correcting for extinction, the source has a bright absolute 
magnitude of 6.99. The absolute magnitude M^ggg and Lx of CX5 place it among the more 
luminous X-ray CVs known in clusters. 



4.5.2. Source CXI 

We discovered the possible counterpart for source CX7 (as well as the counterparts for 
sources CX9 and CX12) by blinking the drizzled U and V images and looking for obvious 
blue sources near the error circles of the X-ray sources. CX7's optical counterpart is clearly 
blue (f/330 — V555=-0.58), although accurate photometry is difficult due to its being located 
in one of the densest parts of the cluster. In order to construct a light curve for this nearly 
blended source, we drizzled the image sets from 2003 into three images, which clearly 
require hours-timescale variability (change of O.Smags) (The star in the 2006 epoch was too 
low in signal in all the images for this procedure to work.) Plotting these data and their 
associated errors shows the star's luminosity is not consistent with a constant source (see 
Fig. [H]). The absolute magnitude of the counterpart in 2006 is My^^g =6.89. In view of 
this object's optical variability and location in the CMD, and the X-ray source properties 



( §3.2.2p . we conclude the optical object is a bright CV and the optical counterpart to source 
CX7. 



4.5.3. Source CX9 



The optical counterpart for this source is very blue (f/330 — V555 = -0.60). It appears 
very bright (Mv555=3.79) in the 2003 V data, but is much dimmer (My^gg =6.065) in 2006. 
It shows variability in the V555 of 0.2mag over the observation in 2003, and 0.6 mag in 
2006 (Fig. [T4|) . From the curve for 2003, we find a lower limit for the period of variation 
of about 8h. Of our counterparts, this is the most convincing case for ellipsoidal variation 
with its small magnitude change, indicating a possible 20-h period. Dramatically, between 
the 2003 epoch and the 2006 epoch, the optical object shows a difference in magnitude 
of more than 2 mags. This is consis tent with dwarf nova eruptions observed in globular 



clusters (e.g. 



Shara et al 



1996 . 



20051). As this object also has similar optical properties to 



the counterpart for source CX5, we use the same argument to identify it as a CV. 



4.5.4. Source CXI 2 

This counterpart to CX12 is very blue (f/sso — V555=-0.49). The counterpart is very 
bright with Mv^^^=Q.Q8 in 2006. Over the 2003 observations it shows variations of 0.2 mag 
and in the 2006 epoch it varies by ~ 0.5 mag (Fig. UM . In light of CX12's CV-like X-ray 
properties, the nature of its optical variation, and its position in the CMD, we conclude 
that it is consistent with a bright CV. 
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5. Discussion 

We find a rich population of X-ray sources in NGC 6388. Our analysis of the density 
profile for the X-ray sources shows them to be more centrally concentrated than the 
brightest stars in the cluster. The inferred characteristic mass of 1.36 Mq for the X-ray 
sources suggests that many of them are likely to be compact binary stars with degenerate 
remnant primaries. Many of these were likely produced in close interactions in the cluster's 
large, dense core. 

We examine some of the characteristics of the different populations of the optical 
sources. The extensive population of blue stragglers was shown to be more massive than 
the horizontal branch stars. The distinctive extreme blue horizontal branch stars were 
shown to be no different in distribution than the other horizontal branch stars, confirming 
the result that dynamical interactions are not a likely cause of the large color distribution 
in the cluster's horizontal branch. 

We identify five sources as bright CVs based on their X-ray luminosities and colors, 
together with the detection of optical counterparts for four of them, and the detection of 
X-ray variability for the remaining one. These counterparts are consistent with bright CVs 
in optical variation, color, and luminosity. In fact, we find at least one of the sources (CX9) 



to exhibit variation on par with that exhibited by dwarf nova outbu rsts. That wou 



Servillat et al. 



d make 



20111). 



it only the 15th such object found in globular clusters to date (see 
The only CV candidate without an identified optical counterpart, CX4, was identified as 
a likely CV based on its X-ray properties; an unusually hard spectrum, well-fit by partial 
covering absorption models, luminosity, and variability. We find five sources that have 
X-ray colors, luminosities, and spectra consistent with what we expect for qLMXBS. These 
sources are not found to vary, which is consistent with expectation for sources without a 
strong power- law component. 
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The CVs found here are some of the brightest (in X-rays and optical hght) yet known 
in globular clusters, suggesting relatively high mass transfer rates. These are likely the 
tip of a much larger CV population. The limit of our detection of CVs was V555 < 23.1, 
for which we detected 4 bright CVs as optical counterparts to X-ray sources. We can 
extrapolate the possible size of the CV population in NGC 6388 by comparison with other 
clusters wi th bright CVs and a de tected CV population. 47 Tuc was found to have at least 



22 CVs by 



Edmonds et al. 



(l2003al ) looking to 8 magnitudes fainter than the MSTO. After 
adjusting for distance and reddening, there are 4 CVs in 47 Tuc with magnitudes similar to 
those in NGC 6388. We can infer that the population of CVs in NGC 6388 is likely as large 
as that of 47 Tuc. This is not surprising given the similarities between the two clusters and 
that they both have some of the highest encounter rates of any cluster (see below). 

The large number of sources detected in our data is consistent with predictions based 
on the encounter rate. We can estimate the encounte r rate using t he approximation given 



m 



Verbunt 



(l2003l ). r oc pirl/vQ. We use values from 



HarrisI (119961) (2010 edition) for po, 



Tc, and vq for this app roximation. 



the clusters studied in 



Pooley et al. 



n Table [71 we give our calculation for this value for 



( I2OO3I ) and for NGC 6388 and compare these to the 
number of detected X-ray sources with Lx > 4 x 10^° erg s~^ in the 0.5-6 keV range. We 
have norm alized our values of F to the mean of the values given for the non-core collapsed 



clusters in 



Pooley et al 



(I2OO3I ). We find them to be roughly consistent, with major 
differences likely due to our use of the updated Harris catalog. NGC 6388 's approximate 
F is seen to be significantly larger than that of all but one of the o ther clusters. In th e 



table we also compare the number of detected sources reported by iPooley et al.l (120031 ) 



with the number detected in more recent studies (as noted in the table comments). The 
larger numbers detected in more recent studies may be attributed to improved reduction 
techniques that allow for a better separation of sources in crowded regions. Ultimately, the 
only two clusters to show a significant difference are u Cen and 47 Tuc, which are most 
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affected by ttie clianges in tlieir values for F, a result of using updated values for the cluster 
parameters. 

In Fig. [151 we plot the number of detected sources, using the newest values when 
possible and correcting for the predicted background contamination, versus F. We note 
that detections for NGC 6388 and NGC 6093 (M80) are not complete to Lx > 4 x 10^° 
erg s~^ (see §2.ip . so their detected source numbers are lower limits (indicated with 



arrows). We fit the relationship of F and N with a po wer law and 



ind it can be described 



Pooley et al. 



( 120031 ). where a similar 



as oc r''-^^='=''-°^. This corresponds with the result of 
trend was found with their calculated values of F. In light of this trend, the detection 
of X-ray sources in NGC 6388 is a confirmation of the prediction that its large, dense 
core produces a large number of close stellar interactions, which in turn accelerate the 
production of close binary systems. 
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Fig. I.- 



Chandra ACIS-S X-ray image of NGC 6388 (0.3-6 keV), wit h the c o re ari d half- 



Harrisl (11996h . 2010 



mass radii indicated (7.2 arcsec and 0.67 arcminutes respectively, from 
version), and the extraction regions for the sources within the half- mass radius of NGC 6388 
marked. 
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Fig. 2. — Representative color X-ray image of the central regions of NGC 6388 (45" x 
45"). North is up and east is to the left. Images in 0.3-1.2 keV, 1.2-2 keV, and 2-6 keV 
were overbinned by a factor of 2, smoothed (using CIAO tool csmooth) using a 0.5" scale, 
and combined using the CIAO tool dmimg2jpg. Likely qLMXBs can be clearly identified by 
their yellow or red colors while the unusually blue color of CX4 stands out, at bottom right. 
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Fig. 3. — Positions of bore-sight corrected X-ray sources overlaid on the t/330 image of 
NGC 6388 from 2006. North is up and east is to the left. Error circles are l-cr errors from 
determination of X-ray source positions; no additional error is included. The large gray circle 
indicates the core radius of 7.2 arcsec. 
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Fig. 4.— X-ray CMD of NGC 6388, plotting X-ray luminosity (0.5-6 keV, corrected for 
extinction of Nh = 2.04 x 10^^ cm~^) against an X-ray color for sources within the half-mass 
radius of NGC 6388. X-ray luminosities are computed using a 6 keV MEKAL spectrum, 
or an average of a MEKAL spectrum and a 10^ K hydrogen-atmosphere NS spectrum, to 
convert photon fluxes into energy fluxes for each of 7 narrow bands, and summed. The 
open squares (colored red in the online version) indicate likely qLMXBs, and open triangles 
(colored blue in the online version) indicate CV candidates with optical counterparts (5,7,9, 
and 12) or identified from X-ray properties (4). 
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Fig. 5. — X-ray spectra (top, data; below, ratio of data/model) of six of the brightest X-ray 
sources in NGC 6388 (a: CXI; b: CX2; c: CX3; d: CX4; e: CX6; f: CX5). We show here 
binned (>10 counts/bin) spectra, though our reported spectral fits use the unbinned data 
and the C-statistic. Sources CXI, CX2, and CX3, which we argue are quiescent LMXBs, 
are plotted with absorbed NSATMOS-F power-law fits. Sources CX4, CX5, and CX6 are 
plotted with absorbed MEKAL fits. Note that CX4 shows an excess at low energies, arguing 
for a partial covering model. 
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Fig. 6. — X-ray light curves for 2 sources, where one source (CX4) shows variation and the 
other (CX2) does not. CX4 is found to vary using the glvary package in CIAO with odds 
ratio 10~^-^ and probabihty 0.043. 
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Fig. 7. — Comparison of the radial distribution of stars brighter than V=20 (the nearly 
smooth hne) with that of all the X-ray sources. The best-fit generahzed King model for the 
X-ray sources is shown by the dotted curve (colored red in the online version). 
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Fig. 8. — Comparison of the radial distribution for non-qLMXBs X-ray sources (the more 
complete distribution curve) vs. qLMXBs identified in this paper. For each distribution, we 
plot the best-fit generalized King model: dotted for the non-qLMXBs and dashed for the 
qLMXBs (colored red and green, respectively, in the online version). 
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Fig. 9. — Comparison of the radial distribution for X-rays sources with counts > 40 vs. those 
with counts < 40. For each distribution, we plot the best-fit generahzed King model: dotted 
for the faint sources and dashed for the bright sources (colored red and green, respectively, 
in the online version). 
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Fig. 10. — CMD for NGC 6388 with 4 optical counterparts labeled by their CX designations 
as used in this paper. The four stellar populations we compare in table [5] and fit in table 
[6] are outlined by the limits used for those populations: RHB (short dashed, colored red in 
the online version); BHB (long dashed, colored blue in the online version); EBHB (dotted, 
colored cyan in the online version); and BS sequence (solid, colored magenta in the online 
version) . 
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Fig. 11. — Finding chart for the center of NGC 6388 in the V555 filter. Er ror circles are 



-a errors on t 



le positions as reported by the authors of A: this resea rch, B 



(2007), and C: 



Liitzgendorf et al. 



( l201lh . For the center labeled D from 



Lanzoni et al. 



Noyola fc Gebhardt 



(120061 ). no error was given by the authors. 
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r (arcsec) 

Fig. 12. — The radial distribution of three different stellar populations. The stars brighter 
than V=20 are depicted by the solid curve, the RHB stars by the short dashed curve and 
the blue stragglers by the long dashed curve. For the BS, we plot the best-fit generalized 
King model (dotted, colored red in the onhne version). Note how the RHB distribution is 
indistinguishable from that of the bright stars, while the BS distribution stands out as being 
more centrally concentrated. 
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Fig. 13. — The finding charts for the four identified optical counterparts. We show the 
F330W from 2006 on the top row, F555W image from 2006 on the middle row, and F555W 
image from 2003 on the bottom row. The circles represent two times the 68% confidence 
regions for the positions of the X-ray sources. The counterparts are noted by the tick marks. 
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Fig. 14. — Light curves with 1-cr photometric error for the counterparts CX5, CX7, CX9, 
and CXI 2 for 2003 and 2006, except in the case of CX7, which was too faint in 2006 for a 
hght curve. 
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Fig. 15. — Plot of detected number of sources vs. the approximate values of F as given in 
Table [71 The source numbers are corrected for background contamination. The dashed line 
represents a least-squares power-law fit. 
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Table 1. HST Data Used in this Study. 



Dataset 


Date 


Filter 


Exposure Time 


GO-9835 


2003 Oct 30 


F555W 


48 X 155 s 




zuuo wet oU 


TT'SI AWT 
r oi4 VV 


lU X ouo s 


GO-9835 


2003 Oct 30 


F814W 


5 X 25 s 


GO-9835 


2003 Oct 30 


F814W 


2 X 469 s 


GO- 10474 


2006 Apr 04/05 


F555W 


48 X 155 s 


GO- 10474 


2006 Apr 04/05 


F814W 


8 X 501 s 


GO- 10474 


2006 Apr 04/05 


F814W 


4 X 25 s 


GO- 10474 


2006 Apr 04/05 


F814W 


4 X 508 s 


GO-10350 


2006 Apr 07 


F555W 


3 X 155 s 


GO-10350 


2006 Apr 07 


F330W 


2 X 1266 s 


GO-10350 


2006 Apr 07 


F330W 


4 X 1314 s 



Note. — All images were captured with the High Reso- 
lution Channel on the Advanced Camera for Surveys. 



Table 2. ST Photometry for Four Optical Counterparts 



star RA (J2000) Dec (J2000) C/330 (2006) ^555 (2003) ^555 (2006) 



5 17'*36"*15.747* -44°44'15.06" 21.99±0.02 22.48±0.02 23.15±0.02 

7 17'*36'"17.187" -44°44'07.57" 22.47±0.08 23.97±0.12 23.05 ±0.31 

9 17''36"*17.231^ -44°44'10.38" 21.63±0.02 19.95±0.01 22.23±0.06 

12 17''36"'17.142^ -44°44'01.79" 22.35±0.04 22.07±0.02 22.84±0.08 



Table 3. X-ray Sources in NGC 6388 





Source 




RA 


Dec 


Distance 


Counts 


Lx J Gr J 


^S S 


Notes 


(CX) 


(CXOGLB J) 


(HH:MM:SS Err) 


(DD:MM:SS Err) 




(0.5-1.5) 


(1.5-6) 


(0.5-6) 


(0.5-2.5) 




1 


173617.6-444416 


17:36: 


17.686 0.002 


-44:44:16.72 0.02 


9.86 


396.0l2i;3 


129.5li?;3 


876.5tll:^ 


842.9ttli 


q 


2 


173616.9-444409 


17:36: 


16.937 0.002 


-44:44:09.85 0.02 


3.42 


268.5li^i 


oo o + ll.l 


663.7tf,:0 


650.2t^^:^ 


q 


3 


173617.3-444408 


17:36: 


17.329 0.001 


-44:44:08.23 0.02 


1.16 


258.81}^;^ 


68.9t^°2' 


796.8ttfl 


700.6+1^:^ 


q 


4 


173616.6-444423 


17:36: 


16.633 0.003 


-44:44:23.48 0.03 


16.50 


34.6t^i 




6Q0.0lTr,l 






5 


173615.7-444415 


17:36: 


15.751 0.003 


-44:44:15.02 0.03 


17.06 




95.21^^ 


444.4t^2-4 


9-1 q ^+24.9 


d 

-J 

1 
1 

c 


6 
7 


173617.3-444406 
173617.1-444407 


17:36: 
17:36: 


17.317 0.002 
17.195 0.003 


-44:44:06.95 0.02 
-44:44:07.59 0.03 


1.67 
0.68 


59.9t?:2 
45.8t|? 




328.6t3g;g 


315.3^3Q g 

173.91^2^:^ 


8 


173618.1-444359 


17:36: 


18.188 0.003 


-44:43:59.47 0.03 


13.53 


97.7t^^/ 


W.9tli 


240.6t^9;2 


240.6t^9;2 


q 


9 


173617.2-444410 


17:36: 


17.245 0.003 


-44:44:10.31 0.04 


2.10 




61.7ti07=^ 


270.8t«;^ 


135.8t?fo 


c 


10 


173617.5-444357 


17:36: 


17.520 0.003 


-44:43:57.13 0.04 


11.54 


86.6t^°3« 




218.1 + n!-? 

— ZO. 1 


214.5+22-q 


q 


11 


173617.0-444403 


17:36: 


17.017 0.003 


-44:44:03.100 0.04 


5.62 


41 7+S-O 

' -6.4 


4i.8t«:li 




128.7ti?;? 


0? 


12 


173617.1-444401 


17:36: 


17.158 0.004 


-44:44:01.76 0.04 


6.48 


33.7t^;? 


41.8t-"3 


264.9t^4-i 




c 


13 


173617.0-444412 


17:36: 


17.088 0.003 


-44:44:12.57 0.03 


4.58 




27.9t^:« 


199.5t^i 


109.5^?^:^ 




14 


173616.8-444412 


17:36: 


16.883 0.004 


-44:44:12.78 0.05 


5.81 




19.8tli 


1+40.1 
iZ0.1_2o 8 


83.7l20;5 




15 


173617.3-444353 


17:36: 


17.346 0.005 


-44:43:53.65 0.06 


14.62 




26.7t^;^ 


1 4+45.4 


rq 7+17.0 




16 
17 


173617.1-444413 
173618.0-444404 


17:36: 
17:36: 


17.158 0.004 
18.054 0.006 


-44:44:13.07 0.04 
-44:44:04.28 0.06 


4.93 
9.73 


14.9^^j 


17 q+6.9 
21.8tl;^ 


144.0t^3.4 
95.9tll^r 


66.21-4 
53.5111;? 





Table 3 — Continued 



Source RA Dec Distance Counts Lx, ergs s ^ Notes 

(CX) (CXOGLB J) (HH:MM:SS Err) (DD:MM:SS Err) (") (0.5-1.5) (1.5-6) (0.5-6) (0.5-2.5) 



18 


173616 


9- 


444411 


17 


36 


16 


969 





005 


-44 


44 


11 


79 





05 


4.43 


q fl+4.8 


18.911? 


100 

lUU.i_2o.9 


40 -|-hl6.8 


19 


173617 


4- 


444403 


17 


36 


17 


405 





007 


-44 


44 


03 


23 





07 


5.36 




14.811^ 


74 


n-i-26.1 

"-15.6 




20 


173617 


5- 


-444406 


17 


36 


17 


558 





007 


-44 


44 


06 


64 





08 


3.94 




9.8^^;^ 


62 


1 -1-32.9 
^-15.0 




21 


173617 


8- 


■444411 


17 


36 


17 


849 





007 


-44 


44 


11 


14 





08 


7.32 




7.8t^;^ 


46 


"-10. 6 


29.7^7^7^ 


22 


173618 


5- 


■444415 


17 


36 


18 


555 





008 


-44 


44 


15 


57 





08 


16.05 


11-61^3 


7.8t^;^ 


39 


q-l-19.2 
"-•-9.3 


26.6^1252 


23 


173617 


3- 


■444357 


17 


36 


17 


322 





007 


-44 


43 


57 


14 





08 


11.12 


5.6t^i 


i3.7t«:^ 


56 


c-l-32.4 
"-14.3 




24 


173617 


8- 


■444344 


17 


36 


17 


863 





008 


-44 


43 


44 


67 





08 


24.52 




I0.8t^;; 


60 


r-l-32.2 
°-16.5 


19.2+_f,,^ 


25 


173613 


9- 


■444415 


17 


36 


13 


911 





009 


-44 


44 


15 


08 





11 


35.91 




7.8t|-J 


51 


q-l-31.3 
''-IS.O 




26 


173617 


4- 


■444410 


17 


36 


17 


440 





009 


-44 


44 


10 


17 





10 


3.06 


8.5t^-i 


4.8t|;S 


27 


1-1-13.8 
■•■-7.6 


2QAtf,i 


27 


173616 


9- 


■444419 


17 


36 


16 


950 





009 


-44 


44 


19 


98 





10 


12.12 




5.8t|-i 


38 


y-l-27.6 

'-12.3 




28 
29 


173617 
173617 


0- 
1- 


■444359 
■444408 


17 
17 


36 
36 


17 
17 


083 
142 






008 
007 


-44 
-44 


43 
44 


59 
08 


86 
70 






09 
08 


8.53 
1.02 


8.8t^-i 

Q 0-I-4.6 
°-°-2.8 


3.9tf:? 
2.9tt-j 


35 
46 


c-l-32.6 
"-12.1 

1-1-35.6 
■■■-14.6 


97 i'-|-14.2 
^'•^-8.9 


30 


173617 


3- 


■444438 


17 


36 


17 


300 





010 


-44 


44 


38 


00 





11 


29.81 






38 


r-l-28.0 
"-•-12.0 


6.51^2:? 


31 


173618 


2- 


-444403 


17 


36 


18 


210 





010 


-44 


44 


03 


20 





11 


11.69 


1.8t?-| 


8.8^2:? 


42 


q-l-32.7 
■^-15.1 


-1 q 0-1-9.5 

^■^•^-5.0 


32 


173618 


4 


444408 


17 


36 


18 


465 





Oil 


-44 


44 


08 


38 





12 


13.28 


2.7t?;^ 


6.9t|^ 


41 


c+33.8 
^-13.9 




33 
34 


173617 
173617 


2- 
2- 


444415 
-444416 


17 
17 


36 
36 


17 
17 


247 
235 






009 
009 


-44 
-44 


44 
44 


15 
16 


09 
86 






10 
10 


6.89 
8.66 


3.8t?;^ 

3-8t?i 


3.9tt° 


21 
23 


c + 16.4 
°-6.9 

r-l-22.1 

^-8.0 


17.4^13^3 



Table 3 — Continued 



Source RA Dec Distance Counts Lx, ergs s Notes 

(CX) (CXOGLB J) (HH:MM:SS Err) (DD:MM:SS Err) (") (0.5-1.5) (1.5-6) (0.5-6) (0.5-2.5) 



35 


173617.3-444421 


17:36 


17.323 0.011 


-44:44:21.27 


36 


173617.0-444351 


17:36 


17.091 0.012 


-44:43:51.02 


37 


173617.6-444408 


17:36 


17.676 0.012 


-44:44:08.84 


38 


173616.2-444415 


17:36 


16.233 0.012 


-44:44:15.43 


39 


173617.6-444404 


17:36 


17.661 0.013 


-44:44:04.89 


40 


173618.9-444358 


17:36 


18.959 0.013 


-44:43:58.08 


41 


173617.0-444405 


17:36 


17.069 0.010 


-44:44:05.94 


42 


173616.8-444406 


17:36 


16.842 0.010 


-44:44:06.59 


43 


173619.0-444341 


17:36 


19.002 0.013 


-44:43:41.79 


44 


173618.5-444412 


17:36 


18.521 0.013 


-44:44:12.24 


45 


173617.5-444402 


17:36 


17.592 0.014 


-44:44:02.16 


46 


173617.5-444413 


17:36 


17.550 0.014 


-44:44:13.88 


47 


173616.5-444410 


17:36 


16.593 0.013 


-44:44:10.63 


48 


173614.7-444430 


17:36 


14.769 0.016 


-44:44:30.79 


49 


173618.1-444419 


17:36 


18.152 0.015 


-44:44:19.24 


50 


173616.0-444408 


17:36 


16.036 0.015 


-44:44:08.14 


51 


173616.6-444421 


17:36 


16.689 0.013 


-44:44:21.94 



12 


13.11 




"•'J-o.o 




I3.5tl;g 


13 


17.24 


6.6t^i 




18.6t-^ 


18.6t-/ 


13 


4.91 


c £+4.1 
^•^-2.2 




23.0l9^4^ 




14 


12.75 




Q 0+4.9 






14 


5.76 




6.8^1^ 




9.51^6 


14 


21.13 




2.9tti 


20.8t^<^6^ 


6.7tli 


11 


2.78 


6.7t^J 






16.5t^°4^ 


12 


4.33 


"•'-2.4 






17.5t^^,° 


13 


32.55 


q 7-1-3.7 


2.8tti 




7 9-1-8.0 
'•^-3.0 


14 


14.45 




3.9tf:0 




9.8t^V 


15 


7.39 


3.6t?:^ 






10.911^72 


15 


6.53 


1.8t?J 


3.8t?:° 


oq -1 -1-26.8 
■^'^■^-9.6 


6.11^:1 


14 


7.09 


4.6t|:^ 






8.4t«:° 


18 


34.53 




3.8l^:° 




"•^-3.0 


17 


14.85 


1.8t?i 


2.9tt:« 






17 


12.61 




2 9+''-^ 


16.9lfg2 




15 


14.85 


0.9t3;° 


2 9+4-'^ 


°-'-4.0 





Table 3 — Continued 



Source RA Dec Distance Counts Lx, ergs s ^ Notes 

(CX) (CXOGLB J) (HH:MM:SS Err) (DD:MM:SS Err) (") (0.5-1.5) (1.5-6) (0.5-6) (0.5-2.5) 



52 


173618.0- 


444412 


17:36 


18.026 


015 


-44:44:12 


70 


53 


173617.3- 


444359 


17:36 


17.363 


015 


-44:43:59 


24 


54 


173616.0- 


444421 


17:36 


16.079 


017 


-44:44:21 


33 


55 


173618.3- 


444434 


17:36 


18.308 


018 


-44:44:34 


22 


56 


173617.0- 


444358 


17:36 


17.074 


015 


-44:43:58 


25 


57 


173618.0- 


444332 


17:36 


18.099 


019 


-44:43:32 


12 


58 


173617.2- 


444341 


17:36 


17.246 


019 


-44:43:41 


74 


59 


173616.7- 


444400 


17:36 


16.752 


019 


-44:44:00 


39 


60 


173616.5- 


444416 


17:36 


16.579 


020 


-44:44:16 


52 


61 


173617.0- 


444329 


17:36 


17.088 


019 


-44:43:29 


88 






16 


9.70 


2 


q+3.5 


1 


0+4.4 









16 


9.10 


2 


q+3.5 

^-1.6 





q+4.4 




5.9t^;? 





19 


17.88 





q+2.9 


2 


0+4.8 
0-1.4 


14.81^^3^ 


0+4.3 
^•'J-1.4 





19 


28.48 


1 


7+3.2 
' -1.1 


1 


q+4.6 
^-1.1 


q + 13.4 







16 


10.08 





q+2.9 

^-0.8 


1 


q+4.6 


9.2+13^5 







21 


37.28 





8+2-9 
^-0.7 


2 


^-1.3 




2.4t|:? 





21 


26.47 


2 


7+3-4 
' -1.4 





n+4.2 
^-0.0 









22 


9.27 


2 


c+3.4 
3-1.4 





n+4.2 









22 


10.75 





7+2.9 
'-0.7 


1 


0+4.6 
^-1.1 


7 6+^2-^ 

' -0-4.3 


q 9 + 6.9 





21 


38.36 





q+2.9 

^-0.7 


1 


p+4.4 


11 7+23.0 
-6.1 


9 4+6.3 



Note. Names, positions, distance from center of NGC 6388, counts in two X-ray energy bands (energies given in keV), 
and estimated X-ray luminosities (in units of 10'^'^ ergs s~^) of X-ray sources associated with NGC 6388. The errors after the 
position represent the la uncertainties in the relative positions of the sources, derived from ACIS_EXTRACT centroiding. 
The counts in each band are the numbers of photons within the source regions of Figure 1. Luminosities are computed from 



the corrected photon fliixes in several narrow bands, see text. Notes indicate short-term variability (V = 99% confidence, 
V? = 95% confidence), years-timescale variability between 2000 and 2003 (Y), and possible identifications (q = qLMXB, q? 
= qLMXB candidate, c = CV with optical counterpart identified, c* = CV identified from X-ray properties, c? = possible 
optical counterpart identification). 



Table 4. Spectral Fits to Bright NGC 6388 Sources 



CX H-atmosphere + Power-law MEKAL Power-law 





kT 


Nh 


PL 
fraction 


Good- 
ness 


Lx 


kT 




Good- 
ness 


Lx 


a 


Nh 


Good- 
ness 


Lx 




(eV) 


(xio^i 

cm-2) 


(%) 


(%) 


ergs s~^) 


(keV) 


(xlO^i 

cm-2) 


(%) 


(xio^^ 

ergs s~^) 




(xlO^l 
cm-2) 


(%) 


( X 10^^ 

ergs s~^) 




1 


I30I2 


4 0+°* 
^■"-0.4 




51 




0.49 


8.9 


100 


62 




8.6i°:l 


8 


130liJ° 


2 


119+1 


=+0.4 
•^■^-0.4 


0+'' 


10 




1 0+0.1 

^••^-0.2 


2.2t°-2 


13.4 


6 7+°-® 

°-'-0.6 


°-"^-0.6 


6.811.1 


2.5 




3 


mtt 


Q 9+0.5 


18111 


40 


1 1 


1.1 


2.2 


95 


7.3 


4 5+°-'^ 
^•^-0.6 


e.olli 


79 


37136 


8 




9 t;+0.5 


Q + 8 


72 




o.8l?,:l 


9 9+0.4 


89 


9 O+0.4 
^■■J-0.3 




o.8t°:3 


24 


c+29 
0_4 


10 


9lti 


9 7+0.7 
^•'-0.5 


n+13 


6 


9 0+0.7 




3.915:'; 


85 


4 0+'-"^ 




7.4ty 


9 


or)+170 
"J '-'-20 


11 




p: c+a.o 
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59 
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4 
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41 


131^ 


80 


7.6 


99 


11.0 


0-9l°o1 




26 


12+_l 
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< 96 


2 ■5+"-^ 
^■^-0.3 
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4 2+1-1 




4.2+ 


24 




9 
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6.4t|« 


4.2i«:^ 


4.011J 
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"^■'-0.8 


9 9+0.6 
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Note. — Spectral fits to sources with more than 80 counts. Errors are 90% confidence for one parameter of interest. The hydrogen column 
density {Njj) to the cluster is 2.2 X lO^^ cm"^^ which was used as a lower bound. Lx in erg s~i for 0.5-10 keV, assuming distance of 10 kpc. 
Hydrogen atmosphere NS fits assume R=10 km, M=1.4 MqNSs, while the power-law photon index was fixed to 1.5. "Goodness" is the fraction 
of simulations with a smaller C-statistic than the given spectral fit; large "goodness" values indicate a poor fit, for which no errors are provided. 
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Table 5. Population distribution comparisons. 



compared pops Prob. 



Optical: 

RHB vs. 20 w/o 0.32 

BS vs. 20 w/o 0.004 

RHB vs. BHB 0.67 

RHB vs. EBHB 0.96 

BHB vs. EBHB 0.91 

X-ray: 

non-qLMXB vs. qLMXB 0.48 

bright vs. faint 0.09 



Note. — The probability that the 
listed populations are drawn from 
the same underlying distribution as 
determined by a K-S test. '20' sig- 
nifies all stars with V555 < 20. 'w/o' 
specifies that the RHB or BS were 
removed from the '20' sample be- 
fore comparing to their respective 
groups. 
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Table 6. Generalized King Model Fit Parameters 



pops 


q 




a 


BS 


1.80 ± 0.27 


4.39 ± 0.48 


-4.39 ± 0.27 


RHB 


0.97 ± 0.10 


7.41 ± 0.93 


-1.92 ± 0.33 


all X-ray with bkd corr 


1.51 ± 0.12 


4.99 ± 0.37 


-3.54 ± 0.37 


qLMXBs 


1.68 ± 0.31 


4.62 ± 0.65 


-4.03 ± 0.92 


non-qLMXBs 


1.40 ± 0.10 


5.33 ± 0.32 


-3.20 ± 0.31 


CVs 


1.57 ± 0.48 


4.85 ± 1.07 


-3.71 ± 1.50 


Bright X-ray (>40 cnts) 


1.82 ± 0.34 


4.35 ± 0.55 


-4.45 ± 1.00 


Faint X-ray 


1.31 ± 0.11 


5.59 ± 0.39 


-2.94 ± 0.33 



Note. — Best-fit generahzed King Model parameters for several 
different populations. Errors are derived from bootstrap statistics. 
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Table 7. Approximate Collision Parameters for a Sample of Clusters 



Cluster 




r 


Sources Detected 


Background 


References 


(1) 


(2) 


(3) 


(4) 




(6) 


(7) 


NGC 6266 


516 


500 


51 




2-3 




47 Tuc 


235 


434 


45 


79 


~ 16 


TT * 1 J_ 1 { C\f\f\r'\ 

Hemke et al. (2005) 


NGC 6626 


148 


186 


26 




2-3 




NGC 6093 


109 


166 


17 


19 


1-2 


Hemke et al. (2003) 


NGC 5904 


40 


69 


16 




4-7 






zo 


zLQ 


28 


60 


1 Q 9zL 


naiitiaiu. ei ai. izuuyi 


NGC 6752 


81 


38 


11 




2-3 




NGC 7099 


23 


18 


7 


11 


1-2 


Lueeer et al. f2007) 


NGC 6121 


8 


13 


6 


6 


1-3 


Bassa et al. r2004) 


NGC 6397 


12 


5.9 


12 


11 


0-1 


BoEfdanov et al. r2010^ 


NGC 6366 


2.1 


2.3 


4 


5 


2-4 


Bassa et al. f2008) 


NGC 6388 


515 






61 


1-2 





Note. — Col. (2): The approximate v alues of F are shown against col. (3): the 



value of r re ported by 



reported by 



Poo 



Pooley et al. 



ev et al 



(I2OO3I ). Col. (4): Numbers of detected sources as 



(120031 ) are compared to col. (5): updated values from the 



literature for a limit of 4 x 10^°erg s ^. The limit for NGC 6388 is ~ 5 x 10^°erg s ^ . 
Col. (6): predicted background counts. Col. (7): sources for the updated results in 
col. (5). 
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